Spectroscopic calculations of the low-lying structure in exotic Os and W isotopes 
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Structural evolution in neutron-rich Os and W isotopes is investigated in terms of the Interacting 
Boson Model (IBM) Hamiltonian determined by (constrained) Hartree-Fock-Bogoliubov (HFB) cal- 
culations with the Gogny-DlS Energy Density Functional (EDF). The interaction strengths of the 
IBM Hamiltonian are produced by mapping the potential energy surface (PES) of the Gogny-EDF 
with quadrupole degrees of freedom onto the corresponding PES of the IBM system. We examine 
the prolate-to-oblate shape/phase transition which is predicted to take place in this region as a 
function of neutron number N within the considered Os and W isotopic chains. The onset of this 
transition is found to be more rapid compared to the neighboring Pt isotopes. The calculations 
also allow prediction of spectroscopic variables (excited state energies and reduced transition prob- 
abilities) which are presented for the neutron-rich 192 - 194 > 196 \y nuclei, for which there is only very 
limited experimental data available to date. 
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I. INTRODUCTION 

Quadrupole collectivity has long been understood as 
one of the most basic, yet prominent, aspects of nuclear 
structure [J Q • Nuclei are quantum many-body systems 
exhibiting collective properties associated with a distinct 
shape of the mean-field, which can be represented by a 
geometrical surface. Quadrupole collectivity can then be 
understood as a quadrupole-sh&ped deformation of the 
nuclear surface whose magnitude depends on the num- 
ber of valence nucleons, and has been shown to exhibit 
remarkable regularities in spectroscopic observables, such 
as the excitation energy of the 2* state and the ratio 
E 4 + / E 2 +. Also evident are stunning shape/phase transi- 
tions at specific nucleon numbers where the collective na- 
ture of the quantal nuclear system can be well described 
as a phase transition between (for example) quadrupole 
vibrational and statically (quadrupole) deformed poten- 
tials HQ. 

The underlying multi-fcrmion dynamics of such nuclei 
however, is so complex that its microscopic understand- 
ing still continues to be a theme of major interest in 
nuclear structure physics research. Mean-field studies, 
based on Skyrme and Gogny Q as well as relativis- 
tic 0,||| Energy Density Functional (EDFs) provide rea- 
sonable descriptions of various nuclear properties such as 
masses, charge radii, mass density distributions, and sur- 
face deformations, over a wide range of neutron and pro- 
ton numbers 0,0]. Such mean-field models, with their 
intrinsic spontaneous symmetry breaking mechanism, are 
highly relevant to understand the microscopy of the nu- 
clear quadrupole deformation and, therefore have been 
used as a starting point for predictions relevant for fu- 



ture nuclear spectroscopic investigations in exotic nuclei 
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Phcnomcnological studies using the Interacting Bo- 
son Model (IBM) [l2j have enjoyed significant success 
in describing the low-lying quadrupole collective states 
of medium-mass and heavy nuclei. The merit of the IBM 
lies in its simplicity, such that, based on group theory, 
the highly complicated multi-fermion dynamics of sur- 
face deformation can be simulated by simple, effective 
bosonic degrees of freedom, which correspond to (collec- 
tive) pairs of valence nucleons [131 ] - In addition to its 
success in reproducing a large amount of experimental 
data on low-lying collective nuclear states in heavy nu- 
clei, the microscopic derivation of the IBM Hamiltonian 
has also been extensively studied (l3l - [l5| . In particular, 
a novel way of deriving the interaction stren gths of an 
IBM Hamiltonian has been proposed recently 16]. This 
method is based on simulating the potential energy sur- 
face (PES) of a given EDF by the corresponding IBM 
PES. The IBM parameters are then derived as functions 
of the nucleon number using the Wavelet analysis method 
[l7l ]. In this way, the universality of the nuclear EDF 
and the simplicity of the IBM can be combined, thereby 
allowing the calculations to predict directly measurable 
spectroscopic observables such as excitation energies and 
electromagnetic transition rates between specific states. 
A number of spectroscopic calculations have been carried 
out using this method for the Ru, Pd, Ba, Xe, Sm isotopic 
chains, as well as theoretical predictions on N > 126 Os- 
W nuclei, using the Skyrme EDF (ITT ] . 

The neutron-rich W, Os and Pt nuclei with A ~ 
190 — 200 exhibit a very challenging structural evolu- 
tion, which has already been extensively studied (l8l - [25| . 
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As originally pointed out in [19( the ratio E A + / E 2 + in 

190 W is anomalously small compared with the one in 
neighbouring isotopes. The most recent experimental 
data on the neutron-rich tungsten chain from 188 > 190 > 192 W 
[2 HH, HH all suggest a change from a well deformed, 
axially symmetric prolate shape for lighter tungsten iso- 
topes, to a more gamma-soft system for 190 W. This tran- 
sition from a prolate to very gamma-soft system for neu- 
tron number N = 116 (i.e., for 190 W) is consistent with 
the recent observation of the second 2 + state in 190 W 
which appears to lie lower than the yrast 4+ in this 
nucleus [21(. The neutron- rich nature of the heavier 
W and Os nuclei make them experimentally challeng- 
ing to study. However, in recent years, there has been 
some progress in their structural investigation following 
multi-nucleon transfer [l|| H3 and isomer and/or 
beta-delayed gamma-ray spectroscopy following projec- 
tile fragmentation reactions pjl, [2l|, [22[ ■ The current ex- 
perimental information is limited to the yrast sequence 
in 190 W [2l|, [H| and the identification of the 2+ state 
in 192 W [2l|. It is interesting to note that the yrast 2+ 
states in the N = 116 isotones 190 W and 192 0s have al- 
most identical energies (~206 keV), as do the N = 118 
isotones 192 W and 194 0s (-218 keV). 

On the theoretical side, mean-field calculations have 
been performed which predict the shapes of these systems 
both with (e.g., [H]) and without (see, e.g., Refs. jHl^. 
and references therein) the assumption of axial symme- 
try in the nuclear mean field. The IBM has also been 
applied to fit the spectral properties of W isotopes in a 
phenomenological way [2!| . More recently, spectroscopic 
calculations have been carried out [3(| to describe the 
structural evolution in Pt isotopes with the Gogny-DlS 
EDF [311 ] . In this paper, we review the current spec- 
troscopy relevant to the prolate-to-oblate shape/phase 
transition in neutron-rich Os and W isotopes. We also 
report the predicted excitation spectra and the transi- 
tion probabilities on the neutron-rich Os and W nuclei. 
The spectroscopic calculations have been carried out in 
terms of the IBM Hamiltonian derived by mapping (con- 
strained) Hartree-Fock-Bogoliubov (HFB) calculations, 
based on the Gogny-DlS EDF, using a similar technique 
as in 1301. 



II. THEORETICAL PROCEDURES 

We begin with the calculation of the PES in terms 
of the (constrained) HFB method using the Gogny-DlS 
EDF. The solution of the HFB equations, leading to the 
set of vacua I^hfb), is based on the equivalence of the 
HFB with a minimization problem which is solved using 
the gradient method [28| . In agreement with the fitting 
protocol of the force, the kinetic energy of the center 
of mass motion is subtracted from the Routhian to be 
minimized, in order to ensure that the center of mass is 
kept at rest. The exchange Coulomb energy is consid- 
ered in the Slater approximation and the contribution of 



the Coulomb interaction to the pairing field is neglected. 
The HFB quasiparticle operators are expanded in a Har- 
monic Oscillator (HO) basis having enough number of 
shells (i.e., N s i ie u — 13 major shells) to grant conver- 
gence for all values of the mass quadrupole operators 
and for all the nuclei studied. The constraint is imposed 
on the average values of the mass quadrupole operators 
Q 20 = \ (2z 2 - x 2 - y 2 ) and Q 22 = & (x 2 - y 2 ) to 
the desired deformation values Q20 = (^hfbI^oI^hfb) 
and Q22 = ($hfb|Q2 2 |*hfb)- In Ref s. [gj the 
Q — 7 energy contour plots with Q = y/ Q 2 Q + Q\ 2 and 
tan 7 = Q22/Q20 have been used to study the (mean- 
field) evolution of the ground state shapes in Pt nuclei. 
Alternatively, one could also consider the (3 — 7 represen- 
tation in which the quadrupole deformation parameter (3 

is written 28] in terms of Q as j3 = \J~^- j^zj 7 where 

(r 2 ) represents the mean squared radius evaluated with 
the corresponding HFB state |$hfb)- The set of con- 
strained HFB calculations provides the Gogny-DlS PES, 
i.e., the total HFB energies £hfb(/3, 7)- 

For the bosonic mapping we use the IBM-2, comprised 
of independent L = + ,2 + proton (s n , d„) and neutron 
(s v , d v ) bosons. The number of proton (neutron) bosons, 
denoted by (n„), is equal to half of the number of 
valence protons (neutrons), assuming the usual magic- 
number shell gaps at Z = 50 and 82, and N = 82 and 126. 
We adopt the standard IBM-2 Hamiltonian flU IetI. f30jj . 
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where h dp = d^ p -d p and Q p = [sj,d p -t-tfj,s p p) +x P [d\d P } { - 2) 
with p = 7r, v. The bosonic PES is represented by the 
expectation value of -ffiBM in the boson coherent state 
;32j , given by 
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where |0) stands for the boson vacuum (i.e., inert core) 
and the coefficients a's are expressed as a P Q = /3 p cos7 P , 
a p ±i — and a p ±2 — -^f3 p smj p . The intrinsic shape 

of the nucleus is then described in terms of the (axially 
symmetric) deformation [3 p and the (triaxial) deforma- 
tion 7p . The IBM PES reads [H,[I3 
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where = j3 v = /3b and 7^ = 7^ = 7b is assumed for 
simplicity [3, E3] • We also assume the proportionality, 
i.e., /3b = Cfl/3, where Cp is a numerical coefficient, and 
7b = 7 [H,[l3]- In this context, the variables /3b and 
7b represent the boson images of the (fermion) defor- 
mation parameters (/3, 7). A point on the HFB PES, 
(/3,7), within an energy range relevant for the considered 
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low-lying quadrupole collective states, is mapped onto 
the corresponding point on the IBM PES, (/3b ,7b)- The 
e, k, x-k,u an d Cp values are fixed for a given nucleus 
by drawing the IBM PES so that the surface topology 
of the corresponding HFB PES is reproduced. This is 
done una mbig uously by means of the recently developed 
procedure [17[ using the Wavelet transform [33j |. 

Note that we compare the total energies -Ehfb(/?,7) 
and Eibm{Ph)- By reproducing the HFB PES as much 
as possible, effects of both vibrational and rotational ki- 
netic energies, similar to those introduced when solving a 
five-dimensional (5D) collective Bohr Hamiltonian (see, 
for example, Refs. 0-HU), should be included in the bo- 
son systems. For large deformation, however, the ro- 
tational response, i.e., the response to cranking, differs 
significantly between nucleon and boson systems, result- 
ing in the deviation of the IBM rotational spectra from 
fermionic ones. This deviation could be corrected by in- 
troducing an additional rotational kinetic- like term, i.e., 
the so-called L ■ L term [34| in the IBM Hamiltonian [35[ . 
This problem does not show up in the present work where 
only the moderately deformed nuclei are concerned, and 
thus we neglect the L ■ L term in the boson Hamiltonian 
of Eq. (P). 

We also note that to what extent the present mapping 
procedure mimics the solution of a 5D Bohr Hamiltonian 
is an interesting open question, which may be partly an- 
swered by looking at how reasonably our results compare 
with the ones of the 5D Hamiltonian and the available 
experimental data. 
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FIG. 1: (Color online) Mapped potential energy surfaces 

(PESs) for !90-196q s and 188-194 W nudei up tQ 3 Me y exd _ 

tation from the energy minimum within the ranges 0° $J 7 ^ 
60° . The PESs are shown in terms of the fermionic deforma- 
tion parameters /3 (= /3b /Cp) and 7 (= 7b). 



III. MAPPED IBM PES AND THE DERIVED 
PARAMETERS 

Figure [1] shows the mapped IBM PESs for i90-i96 Os 
and 188 - 194 W nuclei up to 2 MeV excitation from the en- 
ergy minimum. The corresponding HFB PESs have been 
reported in Fig. 3 of Ref. [Hj]. The PESs for both Os 
and W nuclei show similar tendencies. There are quan- 
titative differences between the Pt and Os-W isotopic 
chains, namely that the topology of the PES changes 
more slowly in the former [301 ] . compared to the latter 
in Fig. rjj An (almost) axially symmetric, oblate min- 
ima is observed in Pt nuclei with TV = 114 ~ 120 and 
shallow triaxiality for N = 110 and 112 [H,|30j]. On the 
other hand, the Os and W isotopes are predicted to have 
the corresponding oblate minima only for N = 118 and 
120, with a more rapid change to axially symmetric pro- 
late deformation for N ^ 114. Indeed, shallow triaxiality 
(i.e., 7-softness) appears only around N = 116 for both 
Os and W nuclei [28|]. The corresponding mapped IBM 
PESs reproduce these trends of the HFB PESs of [H| 
well, whereas the location of the minimum in the IBM 
PES differs from that of the HFB PES of Ref. [H] in 
some nuclei as the present IBM PES of Eq. © does not 
produce a triaxial minimum. The mapped PES for the 
N = 116 isotone, 192 Os is predicted to be very flat along 



the 7-direction. Similarly, the IBM PES for 190 W is also 
very flat, with the global energy minimum corresponding 
to a quadrupole deformation of j3 ~ 0.15 on the oblate 
side. This flatness is the consequence of the x-k and \v 
parameter values, such that their sum is close to zero. 
Comparing Os and W isotopes with the same neutron 
number, the W nuclei are generally steeper in both (3 
and 7 directions than the corresponding Os isotone. A 
similar trend is also observed in the corresponding HFB 
PESs 

Figure [5] shows the evolution of the derived IBM pa- 
rameters for the considered Os and W nuclei as functions 
of the neutron number N. The parameter values for Pt 
nuclei, taken from Ref. 3C}, are also shown for compari- 
son. There are significant differences in quantitative de- 
tails of the derived IBM parameter values between Os-W 
and Pt nuclei. In particular, the values of the parameter 
e in Fig. HJa) for Os and W nuclei are rather small in 
the region away from the shell closure as compared to Pt 
nuclei. In Fig. &h) , the magnitude of the parameter n 
is smaller than the analogous results for the Pt isotopes. 
The behavior of the parameters e and k is reflective of the 
HFB PESs for Os and W nuclei being somewhat steeper 
in the (3 degree of freedom compared to the Pt isotopes, as 
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FIG. 2: (Color online) Derived IBM parameter values for the 
considered Os and W nuclei, represented by solid and dotted 
curves, respectively, as functions of N. Results for Pt isotopes 
taken from Ref. [3y| are also depicted for comparison. 



discussed in Ref. |28j |. The Xn,v parameters in Figs.[2][c) 
and HJd) (as well as their sum) behave similarly to those 
of Pt nuclei. For both 7V=110,112 the sum is almost zero 
in Pt isotopes while it is small for Os and W ones, but 
has a negative sign. This indicates a weak prolate defor- 
mation in the latter as seen in Fig.[TJ In other words, the 
7-soft structure is rather sustained in these Pt isotopes, 
but it is not for the corresponding Os and W isotopes. 
As in Fig. HJe), the scale parameter Cp in the present 
case behaves similarly as for the Pt nuclei with about 
the same order of magnitude. 



IV. CALCULATION OF THE ENERGY 
SPECTRA AND S(E2) VALUES 

Using the derived parameters, we calculate excitation 
spectra and reduced E2 transition probabilities £?(E2). 
The Hamiltonian of Eq. (fTJ is diagonalized by using the 
code NPBOS 0. 

Figure |3l shows ground-state (g.s.) band and the quasi- 
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FIG. 3: (Color online) Low-lying g.s. band and quasi-7- 
bandhead (2+) energies (curves) for 186 - 198 0s and 184 " 196 W 
nuclei. Experimental data (symbols) are taken from Refs. [2ll . 
13, ||. Gogny-DlS EDF is used. 



7-bandhead 2j" (denoted by 2+) energies for Os and W 
isotopes. In general, the calculated results follow the 
experimental trends reasonably well, particularly for 2+ 
energy. What is of interest in Fig. [3] is the behavior of 
the 2+ energy, exhibiting a kink for both 192 0s and 190 W. 
The experimental 2+ energy in 192 Os is lower than the 
4^ one. This is an evidence that the A = 192 nucleus is 
the most 7-unstable one among other Os isotopes. The 
present calculation follows the trend for Os isotopes well, 
and predicts a similar one for W isotopes exhibiting, how- 
ever, more rapid change as a function of N. The location 
of the 2+ state for 196 Os ( 192 W) has not yet been fixed 
experimentally but the present calculations suggest that 
the 4+ state is lower than the 2+ one in both 196 Os and 
192 W. The calculated 2+ energy is generally higher than 
the experimental one, whereas the qualitative feature of 
experimental level is reproduced well. 

Now we turn to the analysis of B(E2) systematics, rele- 
vant to the considered low- lying states. The B(E2) value 
is given by 



B(E2; J -> J') = 



1 



2J + 1 



\{J'\\f^\\J)\ 



(4) 



where J and J' are the angular momenta for the initial 
and final states, respectively. The E2 transition opera- 
tor f (E2 ) is given by f (E2 ) = e^Q v + e v Q v . Here e w 
and e„ stand for the boson effective charges. These ef- 
fective charges should be in principle determined not at 
the mean-field level, but rather by some treatment tak- 
ing into account effects beyond the mean field, such as 
core polarization. This is, however, beyond the scope of 
the current framework and may need to be investigated 
in the future. In what follows, we assume e w = e v , for 
simplicity, and discuss ratios of _B(E2)s rather than their 
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FIG. 4: (Color online) Theoretical (curves) and experimental 
(symbols with error bar) [3^] B(E2) ratios for Os and W iso- 
topes as functions of N. Theoretical results for Pt isotopes 
taken from Ref. [3(J are also depicted as dashed curves, for 
comparison. Gogny-DlS EDF is used. 



absolute values and the quadrupole moments for the cor- 
responding excited states. Note that the B(E2) ratio at 
each dynamical symmetry limit, shown below, means the 
one with infinite boson number [l2T |. 

From Fig. 0|a), we observe that the ratio R\ = 
B(E2; 4+ -> 2+)/B(E2; 2+ -4 Of) does not change much, 
being close to its 0(6) limit of IBM 10/7 (which is also 
the SU(3) limit of R x ). This trend persists for N ^ 118 
where there is currently no available data. The ratio 
R 2 = S(E2;2+ -> 2+)/B{E2;2+ -> 0+), shown in 
Fig. EJb), is of particular interest as one can observe 
a significant difference in its value for the Pt and Os- 
W isotopes. The magnitude of the i?2 ratio is arguably 
the most appropriate and sensitive fingerprint for 7 soft- 
ness jsjj. The i?2 values for both Pt and Os-W are rel- 
atively large and close to the 0(6) limit (= 10/7) for 
N = 114 ~ 118, where the nuclei show notable 7 insta- 
bility. For Pt nuclei, this trend persists even for N ^ 112, 
while smaller values are suggested for Os and W nuclei. 
These differences between the Pt and Os-W chains reflect 
the difference in the topology of the PES. The results for 
Os nuclei follow the experimental trend, which increases 
for N = 110 — 116. The present calculation for Os nu- 
clei suggests the decrease of the i?2 value for N ^ 118, 
which corresponds to a suppression of 7 softness. The ra- 
tio R 3 = B(E2;0 2 f -> 2+)/S(E2;2+ -> 0+) in Fig. Hfc) 
generally has a predicted value which is rather small, be- 
ing close to zero (corresponding to the 0(6) and SU(3) 



limits), as compared to R\ and R2 values. Note that the 
scale of the vertical axis in Fig.Utc) is different from those 
of Figs.lUJa) and@|b). No rapid change with N is seen for 
i?3 as in i?2- Nevertheless, we should note the quantita- 
tive differences between the Pt and the Os-W nuclei. The 
branching ratio R4 = B(E2; 2^ -> 0f)/B(E2; 2% ->■ 2f) 
in Fig. 0Jd) for Os follows the experimental trend for 
N = 110 ~ 116. The decrease of R 4 value from N = 110, 
close to 7/10 (SU(3) limit), toward N — 116, close to zero 
(0(6) and U(5) limit), reflects the corresponding struc- 
tural evolution. The R4 value for the Pt chain is close 
to zero, while for the Os-W chains, there is a signifi- 
cant change at N — 116. For the W nuclei, the ratio 
i?4 increases more rapidly than for the Os chain from 
N = 116 to 112. Earlier phenomenological studies sug- 
gested a similar increase }29j . 

Finally, we present in Fig. [5] the level schemes cor- 
responding to the neutron-rich nuclei 19n ,i92,i94,i96Q s 
and 19 °d92,i94^y taken as representative samples. For 
190 ' 192 Os, for which there are significant experimental 
data, not only the g.s. band but also both the quasi- 
7-bandhead 2+ and the quasi- /3-bandhead 0^ (denoted 
by Ot) energies are reproduced quite well by the cur- 
rent calculations, although the detailed 'in-band' energy 
staggering looks different between the calculated and the 
experimental levels. The calculated B(E2) values for 
190,192q s have been normalized to the experimental [39| 
£?(E2;2+ — ¥ 0+) value. Some algebraic feature is also 
apparent in the calculated results. The At = ±1 rule 
for the E2 decay pattern at the 0(6) limit fl2|, (i.e., 
the dominance of 2\ — > 2\ (Oj — > 2%) over 2\ — > 0+ 
(Oj — > 2^)) in the present calculation also compares well 
with the experimental decay pattern. 

The experimental value of the 2+ energy for 192 Os is 
very close to that of its isotone 190 W (i.e., E «207 keV). 
Also, the excitations energies of the 2+ levels in these 
isotones are also quite similar to each other. The present 
calculations reproduce this overall trend well. In fact, 
the calculated E(2+) = 0.280 (0.278) MeV and 0.286 
(0.274) MeV for 192 0s ( 194 0s) and 190 W ( 192 W) nuclei, 
respectively. For 192 Os and 190 W nuclei, the calculated 
g.s. band energies are rather stretched, and the 2+ en- 
ergies are in good agreement with the respective exper- 
imental data. In the calculated quasi-7 band of Os 
and 192 0s nuclei, one observes a staggering as 2+ (3+ 
4+) (5+ 6+), ... etc.. By contrast, the experimental 
energy spacing shows a more regular pattern. This de- 
viation may be related to the topology of the mapped 
IBM PES in Fig. [TJ which is flat in 7 direction, while the 
correspon ding Gogny-DlS PES exhibits shallow triaxial 
minimum |28[ . In the future, some additional interaction 
term, such as a so-called cubic term [40j . may need to 
be introduced in the boson Hamiltonian to correct the 
deviation for detailed structure of quasi-7 band. 

For 19 4,i96q s nuc i e i ; the predicted 2^ and 4+ energies 
reproduce the experimental ones. The quasi- /?-bandhead 
energy for 194 0s in the present calculation is notably 
larger than the experimental value, which is a conse- 
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FIG. 5: (Color online) Level schemes for ^,^2,194,196 Qg and I90,i92,i94 w nuclei _ The theoretical B(E2) values (in Weisskopf 
units) for 190,192 Os are normalized to the experimental [3j| -B(E2;2+ — ^ 0^) value. Calculated and 4j" energies for 192 Os 
and 190 W are shown (in MeV units) in the parentheses because otherwise these energies look identical. Note that, for 192 W, 
experimental data are shown in the same panel as theoretical ones. Gogny-DlS EDF is used. 



quence of the peculiar topology of the Gogny-HFB PES, 
which exhibit a pronounced oblate minimum with a rela- 
tively small deformation. This results in the larger value 
of the parameter k than the one in the IBM phenomenol- 
ogy [29[ which would give good agreement for the excited 
+ energies. The positions of the 2+ and the 0t energies 
for 196 0s are predicted to lie below and beyond the 6^ 
level, respectively. For the exotic 192 W and 194 W nu- 
clei, the present calculation suggests a quite similar level 
pattern to their respective isotones, 194 0s and 196 0s. 



V. SUMMARY 

To summarize, we have presented the predicted ex- 
citation spectra and B(E2) ratios of exotic Os and W 
isotopes with N = 114 ~ 120. Spectroscopic calcula- 
tions have been carried out in terms of the IBM Hamil- 
tonian constructed by the constrained HFB calculations 
with Gogny-DlS EDF. We have examined the prolate- 
to-oblate shape/phase transition as functions of neutron 
number N in the considered isotopic chains. The exper- 
imental trends of not only g.s.-band energies but also 
the quasi-7-bandhead 2+ energy for Os isotopes is re- 
produced well, suggesting that the N = 116 nucleus is 



the softest in 7. A similar pattern is predicted in W 
isotopes, while the evolution of levels appears to occur 
more rapidly in W than in Os. Interestingly enough, all 
these results reflect to a good extent the results of the 
underlying microscopic Gogny-HFB calculations. Lastly, 
let us comment on the form of the boson Hamiltonian in 
Eq. ([1]). While this form may be rather simple, it deter- 
mines the basic topology of the PES, and is supposed to 
be the most relevant for the description of the low-lying 
structure at the present stage. On the other hand, the 
IBM-2 phenomenology considers additional interaction 
terms as compared to those in Eq. ([1]). Some of these 
terms have a minor effect, but others might affect the 
spectroscopic results in quantitative details as suggested 
in the structure of quasi-7 bands in Fig. [5j It should be 
then very interesting to study in the future which parts 
of a more general boson Hamiltonian are crucial, as well 
as how they affect the spectroscopic properties quantita- 
tively. 
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